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We investigate a chiral property of the domain- wall fermion (DWF) system using the four-dimensional hermitian 
Wilson-Dirac operator Hw- A formula expressing the Ward-Takahashi identity quark mass ms, with eigenvalues of 
this operator is derived, which well explains the dependence of ra^q observed in previous numerical simulations. 
We further discuss the chiral property of DWF in the large volume in terms of the spectra of Hw- 



1. Introduction 

Recently a suitable characterization of chiral 
symmetry on the lattice has appeared in the form 
of the Ginsberg- Wilson relation [|l|. A num- 
ber of numerical analysis have been carried out 
for the domain wall fermion (DWF) ||^, which 
is one of the simplest solutions of this relation 
lll^. A chiral property of DWF is investigated 
numerically in two ways; using either the chi- 
ral symmetry breaking Ward-Takahashi identity 
mass |^,^ or the zero-eigenvalue density p(0) 
of the four-dimensional hermitian Wilson-Dirac 
operator Hw The latter is also related to 

the parity-flavor breaking order parameter |^,|lO| . 
However, conclusions are contradictory whether 
chiral symmetry is realized in DWF. 
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In this article we derive a formula which ex- 
presses m^q in terms of the eigenvalues of Hw-, 
which helps to resolve the mutual inconsistency. 
We further discuss the chiral property of DWF, 
investigating the lattice volume dependence of 



in this formula. 



2. Formula 

The anomalous quark mass m^q is defined 
through a Green function of the chiral symme- 
try breaking term Jsg 1^ of the Ward-Takahashi 
identity: 



(1) 



where P{x) is the pseudo scalar density. 

We try to express ni^q in terms of the eigen- 
values of Hw- For this purpose we expand the 
Green function by inserting 1 = ("-I where 

\n) is the eigenstate of a deformed Hamiltonian 
H defined by 
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Figure 1. J2x I (^^l") P as a function of t. 

To simplify the resulting expression, we no- 
tice that the eigenvalues of H\y or H are clas- 
sified into isolated eigenvalues distributed in the 
lower region and almost continuous ones in the 
upper region. As is seen in Fig. |l| the eigen- 
functions for the isolated eigenvalues are local- 
ized exponentially around a center and higher ex- 
cited states with continuous eigenvalues are al- 
most plane waves, rapidly oscillating before ab- 
solute square is taken. 

We adopt an exponential localization for the 
low modes and a plane wave function for the ex- 
cited modes as an approximation. The formula 
for TOsg then becomes 



TO5g = 
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with V the space-time volume, A an eigenvalue of 
H and h being related to the exponential width 
6 of the localized eigenfunctions as h—{46)'^. 

3. Numerical results 

In Figs. ^ and |^ we compare m^q obtained with 
(^) (lines) with those directly calculated in sim- 
ulations, both for the plaquette and the RG im- 
proved gauge actions, where circles and squares 
are from ||^,^ and triangles are from this work. 

The right hand side of (||) is calculated as fol- 
lows. We evaluate eigenvalues of the Wilson Dirac 
operator Hw using the Lanczos method. The 
eigenvalue A of the deformed Hamiltonian is eval- 
uated from that of Hw perturbatively for the iso- 



Figure 2. m^q as function of N5 for the plaquette 
action. Lines are our formula. RBC is a data 
from [pi. 
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Figure 3. as function of for the RG action. 

lated eigenvalues and directly for the continuous 
ones assuming the free theory relation. The sep- 
aration between the isolated and the continuous 
eigenvalues is made byjnspection of the eigen- 
value distribution, and h is treated as a free pa- 
rameter to fit the data. For example, at /3 = 2.6 
of the RG action, we adopt the number of local- 
ized eigenvectors n/ = 20, and h = 200. Summing 
up contributions from each eigenvalue we obtain 
m^q as a function of N^. Errors are calculated by 
a single elimination jackknife procedure for 100 
configurations. 

As is seen in these figures our formula well ex- 
plains the iVs dependence of the anomalous quark 
mass both in strong and weak coupling regions for 
the two gauge actions. An important problem in 
the previous numerical investigations is whether 
TOsg contains a constant term in the large 
limit. Now we can conclude that the exponen- 
tial dumping still holds even at large but with 
a small decay rate, which reflects the magnitude 
of a few small eigenvalues. This is a reasonable 
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Figure 4. Volume dependence of 7715, for the pla- 
quette action at (3 = 6.0. 
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Figure 5. Volume dependence of m^q for the RG 
action a,t (3 — 2.6. 



behavior for finite lattice volume where exact zero 
eigenvalues are absent. 

A bending behavior of 7715^ for weak coupling 
around TVs ^ 20 can be explained with our formula 
as follows. For small the continuous modes, 
which dominate in number, mainly contribute in 
(H), giving a steep exponential decay. For large 
A^5, on the other hand, only a few of the small iso- 
lated eigenvalues contribute, so that the dumping 
rate of m^q becomes small. 

4. Volume dependence of m^q 

Although an exponential decay of m^q in Nr, 
seems to hold in DWF at finite volume, there 
remains a crucial question whether this remains 
so, and hence chiral symmetry is exactly real- 
ized with DWF, in the infinite volume limit. In- 
deed TOsg decreases with some power of I/N5 if 
the number of small eigenvalues increases linearly 
with the volume V 111]] , in which case chiral sym- 
metry is not realized. 

To examine this question we plot the volume 
dependence of at a weak coupling in Figs. ^ 
and ^, using the formula (|^) with n; and h ob- 
tained in the previous section. Almost no volume 
dependence is observed for V > 12^ at large 
for both actions, and a similar behavior holds also 
in the strong coupling region. 

It is difficult, however, to distinguish a slow 
exponential decay seen in Figs. || and || from a 
power of l/iVs at ^ 100. A distinction be- 
comes possible only for ^ 100, which effec- 
tively corresponds to iVs ^ cx) for the current 
lattice volumes, V ~ 12^ - 24*, so data at much 



larger volumes are needed to establish the stabil- 
ity of an exponential decay against the increase 
of volume. Therefore the current data do not al- 
low definite conclusions on the chiral property of 
DWF in the infinite volume limit. Further inves- 
tigations will be required to answer the question. 
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